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ABSTRACT 

The broad line region (BLR) of luminous active galactic nuclei (AGN) is a prominent observational 
signature of the accretion flow around supermassive black holes, which can be used to measure their 
masses (Mbh) over cosmic history. Due to the < 100 /xas angular size of the BLR, current direct 
constraints on BLR kinematics are limited to those provided by reverberation mapping studies, which 
are most efficiently carried out on low-luminosity L and low-redshift z AGN. We analyze the possibility 
to measure the BLR size and study its kinematic structure using spectroastrometry , whereby one 
measures the spatial position centroid of emission line photons as a function of velocity. We calculate 
the expected spectroastrometric signal of a rotation-dominated BLR for various assumptions about 
the ratio of random to rotational motions, and the radial distribution of the BLR gas. We show that 
for hyper-luminous quasars at z < 2.5, the size of the low-ionization BLR can already be constrained 
with existing telescopes and adaptive optics systems, thus providing a novel method to spatially 
resolve the kinematics of the accretion flow at 10 3 — 10 4 gravitational radii, and measure Mbh at 
the high-L end of the AGN family. With a 30m-class telescope, BLR spectroastrometry should be 
routinely detectable for much fainter quasars out to z ~ 6, and for various emission lines. This will 
enable kinematic Mbh measurements as a function of luminosity and redshift, providing a compelling 
science case for next generation telescopes. 


1. INTRODUCTION 

Quasars are the most luminous compact objects in the 
Universe, with luminosities L reaching up to 10 48 erg s _1 . 
Such immense energy outputs are believed to be the re¬ 
sult of gas accretion onto black holes (BHs) with masses 
of up to Mbh ~ 10 10 M Q (Lynden-Bell 1969), shin¬ 
ing at luminosities reaching their Eddington luminosity 
(1.3 x 10 38 Mbh/Mq erg s _1 ). 

Arguably the most distinct property of quasar spec¬ 
tra is the prominence of broad emission lines, with typ¬ 
ical widths of ~ 3,000 km s _1 (e.g. Vanden Berk et al. 
2001), and maximum widths of 10 — 20, 000 km s _1 (Laor 
2003; Fine et al. 2008; Stern & Laor 2012a). The ioniza¬ 
tion structure of the broad line emitting gas, known as 
the broad line region (BLR), is well described by photo¬ 
ionization by the central continuum source (Davidson & 
Netzer 1979 and citations thereafter). However, the spa¬ 
tial structure and dynamics of the BLR are not well un¬ 
derstood. The large observed velocities suggest that the 
kinematics of the BLR are dominated by the BH gravity, 
and that the BLR resides at ~ 10 3 — 10 4 r g , where r g is 
the gravitational radius. 

Assuming the observed BLR velocities are of the order 
of the local Keplerian velocity, then a measure of the BLR 
distance from the BH, tblr, provides a measure of Mbh- 
Since luminous quasars are observable to high redshift z, 
then estimates of Mbh in quasars can be utilized to track 
the buildup of massive BHs over cosmic time. In low-L 
and low-2 Active Galactic Nuclei (AGN), the family of 
objects of which quasars are the luminous subset, tblr 
can be measured with reverberation mapping (RM). In 
RM, one derives tblr from the time lag between changes 
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in the continuum luminosity and the corresponding vari¬ 
ation in the line luminosity (Blandford & McKee 1982; 
Peterson 1993; Peterson et al. 2004). RM has been widely 
applied in the past two decades, and has yielded three 
main results. First, the response of H/3 to changes in 
continuum luminosity in single objects suggests that the 
BLR resides in a narrow range of r, where the bulk of the 
H/3 emission originates from a dynamical range of 5 — 10 
in r (Maoz et al. 1991; Pancoast et al. 2013). Second, 
the characteristic response-weighted tblr of H/3 scales 
as ~ L 1//2 (Kaspi et al. 2005; Bentz et al. 2009b, 2013), 
with a typical value of 0.05 pc at L45 = 1, where L45 is 
the monochromatic luminosity at 1450A L 1450 ^ in units 
of 10 45 erg s _1 . And third, different broad emission lines 
are emitted from somewhat different r, where generally 
higher ionization lines originate from smaller r than the 
Balmer lines. 

The observed value of tblr, its scaling with L, and 
its small dynamical range, are all aptly explained by 
a combination of two effects which suppress the line 

emission at r <C 0.05 L45 pc and r > 0.1 L45 pc. At 

r > r su b « 0.1 L 45 pc, where r su b is the dust sublimation 
radius, dust grains can survive and hence suppress the 

line emission (Netzer & Laor 1993). At r -C 0.05 L 4 ' 5 pc, 
the equilibrium of BLR gas pressure with radiation pres¬ 
sure suggested by Baskin et al. (2014) implies that the 
gas volume densities are so high (j> 10 11 cm -3 ), that 
emission of the permitted lines is collisionally suppressed. 
The equilibrium with radiation pressure can also explain 
the stratification of lines with ionization (see fig. 5 and 
§4.6 in Baskin et al. 2014). 

Do the aforementioned properties of tblr apply also 
to high-L and high-z quasars? RM of a couple of 
L « 10 47 erg s _1 quasars at z ~ 2 suggests that the 
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answer is yes (Kaspi et al. 2007; Chelouche et al. 2012). 
However, applying RM to luminous quasars is problem¬ 
atic since in high-L quasars telr is expected to be on the 
scale of light-years, with correspondingly long response 
times. Also, since quasar evolution dictates that high-L 
quasars are at high-z, cosmological time dilation makes 
the observed time lags even longer. Furthermore, high 
luminosity quasars are less variable (Angione & Smith 
1972; Vanden Berk et al. 2004) and therefore have a 
weaker RM signal. The combination of a weak signal 
with RM timescales of several years in the observed frame 
makes RM of high-L high-z quasars extremely challeng¬ 
ing, and hence to date RM studies of only a couple of 
such objects have been published. Therefore, an inde¬ 
pendent method to constrain telr (and Mbh) at high-L 
and high-z would be a big advantage. 

Another open question concerns the kinematic struc¬ 
ture of the BLR, which given that it resides at ~ 10 3 — 
10 4 r g , is likely an integral part of the accretion flow. 
Is the BLR in some kind of ordered flow, such as an 
extension of the accretion disk, or are the motions dis¬ 
ordered, such as a population of clouds in random virial 
motion about the BH? There are some indirect observa¬ 
tions which favor an ordered flow, such as the rotation 
of the polarization angle across the Ha spectral profile 
seen in some objects (Smith et al. 2005), and the distinct 
variability characteristics of the red and blue line wings 
which is sometimes observed (Veillcux & Zheng 1991). 
Both the rotation of the polarization angle and the dis¬ 
tinct variability patterns as a function of velocity suggest 
that the red-wing emitting gas and the blue-wing emit¬ 
ting gas are spatially distinct, as expected in an ordered 
flow. Also, an ordered flow is favored due to the lack 
of ‘spikes’ in the broad emission line profile, which are 
expected if the BLR is composed of individual clouds in 
a random velocity field (Arav et al. 1998; Dietrich et al. 
1999; Laor et al. 2006). However, a direct observation 
which confirms the existence of an ordered velocity field 
in the BLR remains elusive. 

In this study, we analyze the possibility to constrain 
the BLR size and velocity field in high-L and high-z 
quasars using spectroastrometry (e.g. Bailey 1998). In 
spectroastrometry, one measures the position centroid 
of photons as a function of photon wavelength, which 
in principle can be pinpointed to an accuracy which is 

~ TVpij higher than the resolution limit of the tele¬ 
scope, where N p h is the number of photons. Extrap¬ 
olation of the 7 ~blr — L relation mentioned above to 
the most luminous quasars implies BLR angular sizes 
of s=s 100 /xas, which is a factor of ~ 10 3 below the resolu¬ 
tion of « 50 mas of diffraction-limited observations on 8 m 
telescopes using adaptive optics at near-infrared wave¬ 
lengths. Therefore, with ~ 10® photons, the 100 /xas- 
scale photocenter offset between the red-wing and blue¬ 
wing photons can in principle be detected. As spec¬ 
troastrometry is most effective in luminous quasars where 
rBLR is large and the photon flux is high, it is comple¬ 
mentary to RM which is most easily applied to low-L 
AGN. Additionally, spectroastrometry measures an r- 
weighted function of the BLR (see below), compared to 
the response-weighted function of the BLR measured by 
RM, hence the two methods give independent constraints 
on the size, kinematics, and distribution of material in 


the BLR. 

Bailey (1998) applied spectroastrometry to seeing lim¬ 
ited observations of binary stars and of the narrow line re¬ 
gion (NLR) in AGN. Spectroastrometry was also applied 
to nuclear gas in the Circinus galaxy (Gnerucci et al. 
2013), to molecular disks in young stellar objects (Pon- 
toppidan et al. 2008, 2011; Joergens et al. 2013), and to 
planetary nebulae (Blanco Cardenas et al. 2014). Using 
Adaptive Optics, Pontoppidan et al. succeeded in over¬ 
coming atmospheric and technical issues and achieved 
photon-limited angular resolutions of 100 — 500 /i as, sug¬ 
gesting that resolving the BLR with spectroastrometry 
is feasible with existing telescopes. 

Conceptually similar to BLR spectroastrometry, Shen 
( 2012 ) discusses the possibility to exploit the astromet¬ 
ric position of a variable broad line emission, via imag¬ 
ing with a narrow band filter with wavelength tuned 
to the broad emission line. As discussed in Trippe 
et al. (2010), measuring the 2D-astrometric position is 
dominated by optical aberrations in the imaging instru¬ 
ment, and current imagers show astrometric biases at the 
400 /ras level, a factor of >4 too large for BLR spectroas¬ 
trometry. Hence, Shen (2012) concluded that detecting 
the astrometric signal of the BLR will be delayed to 
next-generation telescopes. However, in spectroastrome¬ 
try the optical abberation bias is circumvented, since all 
photons practically travel along the same path through 
the optics (see below), and therefore the astrometric sig¬ 
nal can in principle be detected already by existing tele¬ 
scopes. Additionally, since the technique suggested by 
Shen (2012) requires a variable BLR, it has the same ob¬ 
servational resource requirements as RM, and is there¬ 
fore practically limited to focus on the same highly and 
shortly variable low-L low-z objects. Therefore, in this 
paper we limit the discussion to the simpler possibility 
of applying spectroastrometry to the BLR without any 
temporal information. 

This paper is structured as follows. In §2 we calculate 
the BLR angular sizes for the most luminous quasars at 
each z. In §3 we estimate the expected spectroastro- 
metric signal for a rotation-dominated BLR, for differ¬ 
ent assumptions on the radial distribution of the BLR 
gas and on the ratio of random to rotational motion. We 
present several observational considerations in §4, and 
simulate the expected spectroastrometric signal with 8 m 
and 30m-class telescopes in §5. We discuss our results 
in § 6 , and summarize in §7. Throughout the paper, we 
assume a FRW cosmology with H = 0.3, A = 0.7 and 
Ho = 70 km s _1 Mpc -1 . 

2. THE ANGULAR SIZE OF THE BROAD LINE REGION 

Kaspi et al. (2005) derived the relation between L 1450 ^ 
and trm, the distance of the BLR measured by rever¬ 
beration mapping of H/3, at a luminosity range 10 41 < 
L 1450 ^ < 10 46 erg s _1 . To estimate 7"blr in spec¬ 
troastrometry candidates, we extrapolate their relation 
to the higher luminosities considered here L 1450j ^ > 
10 47 erg s _1 : 

r »“ = O5 ( 1 O«erg > S -0 pC ■ (1) 

We can convert eqn. (1) to a cosmology-independent 
relation based on angular size 6 and observed flux. De- 
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Figure 1. The apparent magnitudes and expected BLR angular 
sizes of the brightest AGN at each z, complied from BQS, SDSS, 
the z > 5.7 quasar sample of Banados et al. (2014), and a sample 
of local Seyferts. (Top) For each 0.5-wide bin in z, the values 
of mi (K-corrected to z = 2) of the brightest quasar, the 10 th 
brightest, and the 100 th brightest quasar, are shown. The ranges 
in z where Paa, Ha, and Mg n are observed in NIR bands are 
noted on top. (Bottom) The expected BLR angular size of the 
quasars shown in the top panel, derived from L 1450j ^ by assuming 
the RM-based relation between trm and L 1450j ^. The implied 
cosmology-independent relation between 6 and F u (eqn. 3) is noted. 
The dependence of the maximum 0(7*rm) on z is remarkably weak, 
spanning merely a factor of three over the entire z range. 

fine F v . 145oA as the flux density at rest-frame wavelength 

Arest = 1450A, and defining d l as the luminosity dis¬ 
tance, we have 

L4 5 oA = TTT^;1450A'4^ ■ (2) 

Eq. (2), combined with the relation between angular dis¬ 
tance and luminosity distance rRM/^(^RM) = ^/(l+z) 2 , 
gives 

%Rm) = 51 ( xOmTy ) ^ + ^ ^ 

It is convenient to express eqn. (3) also with the com¬ 
monly used nrii{z = 2), the apparent i-band magnitude 
after K-correcting to z = 2 (i.e. the apparent magnitude 
at A re st = 2500A). Therefore, 

F v . 145oA = 0.74 • 10-0.4(^21+48.6) } (4) 

where the factor of 0.74 follows from a Richards et al. 
(2006) quasar template. Together with eqn. (3), we get 

0(r RM ) = 22 • io( 15 -™i( J = 2 ))/ 5 (l + z f/ 2 ^as . (5) 

We search for the brightest quasars at each z in the 
Bright Quasar Survey (BQS, Schmidt & Green 1983), 
the quasar catalogs of the Sloan Digital Sky Survey 7 th 
and 10 th data releases (SDSS-DR7, Schneider et al. 2010; 
SDSS-DR10, Paris et al. 2014), and the z > 5.7 quasar 
sample compiled by Banados et al. (2014). We sup¬ 
plement this quasar list with local bright Seyferts from 


Bentz et al. (2009a) , who corrected the AGN continuum 
emission for the host galaxy contribution, which may be 
significant in the relatively low luminosity Seyferts. We 
derive m\{z = 2) for SDSS quasars by K-correcting the 
PSF magnitude of the SDSS-band which is closest to 
2500A in rest-frame wavelength. For BQS quasars we de¬ 
rive rrii(z = 2) from the B-magnitude, and for the high -2 
quasars we use the J-band magnitudes, if available, or 
otherwise the y-band magnitude. The luminosity of the 
local Seyferts is based on the luminosity at 5100A listed 
in table 9 of Bentz et al. (2009a). All K-corrections are 
performed assuming a Richards et al. (2006) quasar tem¬ 
plate. We divide the quasars into bins of 0.5 in 2 , and 
sort them by brightness. The brightest 2 , 10 th brightest, 
and 100 th brightest quasar in each 2 -bin are shown in 
Figure 1. For reference, we note the ranges in 2 where 
Paa, Ha, and Mg II are observed in NIR bands. As ex¬ 
pected, nil (2 = 2 ) of the brightest quasars generally de¬ 
creases with 2 . However, due to the factor of (1 + 2) 3//2 
in eqn. (5), the dependence of the maximum 0 (t R m) on 
2 is remarkably weak, spanning merely a factor of three 
(50 — 150 /Lias) over the entire 2 range. 


3. THE SPECTROASTROMETRIC SIGNAL OF THE BLR 


Chen et al. (1989) and Chen & Halpern (1989, here¬ 
after CH89) derived the spectral profile of a broad emis¬ 
sion line originating from a rotating disk, including both 
special and general relativistic effects to first order in 
r s /r. The two main parameters of their model are the 
distribution of line emission as a function of r, and the 
amount of line broadening at each location of the disk, 
due to unordered motions which could be sourced by tur¬ 
bulence or other physical processes. Below, we discuss 
the physical and observational constraints on these two 
parameters. Here, we use the CH89 formulation to de¬ 
rive the expected spectroastrometric signal of a rotating 
BLR as a function of these two parameters. 

We note that the flat disk assumed by CH89 is prob¬ 
ably an oversimplification even for a rotation-dominated 
BLR, since the BLR characteristic height above the disk 
mid-plane may change with r. Therefore, our analysis is 
based on the assumption that as long as the dominant 
motion is rotation, the CH89 formulation should give a 
reasonable approximation of the expected spectroastro¬ 
metric signal. 

Following CH89, we define r and tp’ as the polar coordi¬ 
nates in the disk frame. We assume that the orientation 
on the sky of the spatial direction of the telescope slit is 
at an angle j from the major axis of the projected BLR 
ellipse. These definitions are pictured in the top panel 
of Figure 2, which depicts an idealized BLR, which orig¬ 
inates from a single r = tblr- 


We define $*(r, tp')drdtp' 


— 2 —1 
cm s 


(km s 1 


as 


the observed flux density of photons per unit velocity 
which originate from the disk coordinate Simi¬ 

larly, the photocenter position of these photons along the 
spatial direction of the slit is defined as S*(r,<p') [/ras]. 
The value of S* (r, tp') can be calculated from the impact 
parameter of a photon at infinity 6 , which equals (eqn. 6 


2 None of these quasars appears in the list of lensed quasars 
compiled by Veron-Cetty & Veron (2010). 
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Figure 2. (Top) An illustration of a spectroastrometric observa¬ 
tion of a simplified BLR, which originates from a single r = tblr 
and has negligible local line broadening. The ellipse depicts the 
BLR ring, with coordinate y/, projected onto the plane of the sky. 
Color denotes line-of-sight velocity. The physical size of the ellipse 
axes, up to relativistic effects, are tblr and ^blr cos L The angle 
between the spatial direction of the slit and the major axis of the 
projected ring is defined as j. (Middle) The implied emission¬ 
line profile of the rotating ring on top of the underlying quasar 
continuum, for an assumed u ro t = 10 000 km s —1 , sinz = 0.5, and 
EW = 26 000 km s —1 . Normalization assumes an Ho; line of a 
L = 10 48 erg s —1 quasar at z = 2. (Bottom) The implied photo¬ 
center offsets of the BLR, assuming no contamination by non-BLR 
sources, for the ^(^elr) and j noted in the panel. The centroid of 
the reddest line photons are offset from the centroid of the bluest 
line photons by 150 p,as. 


in CH89) 

b = r (l — sin 2 i cos 2 y /) 1 ^ 2 (l + 0(— )) , ( 6 ) 

where i is the inclination of the disk normal to the line 
of sight. The first term in eqn. ( 6 ) is the Newtonian 
projection of (r, tp') on the sky, and the second term is 
a correction due to light bending by the central mass. 
Projecting b on the slit spatial direction yields, with some 
trigonometry, 

S*(r, ip') = r (sin j cosicosip' + cos j siny/) (l + 0(— )) . 

( 7 ) 

The observable photocenter position at velocity v, S v , 
is the photon-flux-average of S* (r, ip'): 

o _ If S Z( r ’P')®Z( r ’ < P') d< P' dr /o\ 

S v -- - - , (8) 

where is the observed photon flux density: 

$v= JJ $*(r, y/)dy/dr . (9) 


Now, since <fr*(r, y/) = 4>*(r, n — ip') (1 + 0(r g /r)) (see 
CH89), then by summing the contribution from ip' and 
7 r — ip' to the integral in eqn. ( 8 ), the term in eqn. (7) 
with cosy/ cancels out to zeroth order in r g /r, and we 
get: 

\ K(r,y/)<I>*(r,y/) + S*(r,n - y>')$*(r, 7 r - ip')] = 

rcosjsiny/<I?*(r, ip') (l + 0 ( —)) ( 10 ) 
Hence, eqn. ( 8 ) simplifies to 


C • f rdr f dip' sin y/$*(r, <p') (l + O(-jr)) 

J v — COSJ • -—- 

The light-bending corrections to S v are of the order 
of 0(r g /r blr) 10~ 3 , and are likely weaker than the 
inaccuracy of the model due to the assumption that the 
BLR is a flat rotating disk, so it is not clear whether 
these corrections are observable 3 . 

The observed S v will be diluted by the quasar contin¬ 
uum, which at optical and ultraviolet frequencies orig¬ 
inates from the accretion disk. The continuum emis¬ 
sion will have a negligible spectroastrometric signal, and 
therefore will dilute S v of the broad line by a factor of 
$„/(<!>„ + < h™ nt ), where 4>™ nt is the photon flux density 
of the continuum 4 . Therefore, the final form of eqn. (11) 
is 


0 . / rdr J dip' sin y/$*(r, ip') (l + O(qr)) 

£> v - COSJ ■ ^ , $cont { LZ ! 




3.1. A Simplified BLR 

In order to explore the dependence of S v on the disk 
properties, we begin by calculating S v for a simplified 

3 Note that relativistic corrections of orders 0((ry/Y blr) 1 //2 ) 
and 0(r g /rBL,R) to the hne P ro fi^ e have been observed (CH89; 
Tremaine et al. 2014). 

4 Equivalent to eqn. 4 in Pontoppidan et al. (2008). 
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BLR which originates from a single thin ring r = tblr 
and has negligible local line-broadening, i.e. a <C v rot , 
where a is the characteristic velocity of the local line¬ 
broadening mechanism, and v ro t is the rotational veloc¬ 
ity. Under these conditions, v rot = (GMbh/lblr) 1 ^ 2 - 
We relax both assumptions below. In the figures, we use 
the fully relativistic expression for <h*(r, <p') from CH89 5 , 
and eqn. (12) to calculate S v . However, since the rela¬ 
tivistic corrections to $*(r, <p') have a weak effect on the 
derived S v , for increased readability we repeat here the 
Newtonian expressions for d>*. In this simplified BLR, 
<f>* is equal to 

3>*(r, ip') oc 6{r — r B LR)6(sin<//-— T ) (13) 

^rot sin % 


The implied is shown in the second panel of Fig. 2, 
assuming v ro t — 10 000km s -1 , i = 30°, and a total 
broad line photon flux of 4 x 10 6 photons hr -1 m -2 . This 
broad line flux is appropriate for the broad Ha line of a 
L = 10 48 erg s” 1 quasar at z = 2, assuming an equivalent 
width (EW) of 570A compared to the AGN continuum 
(Stern & Laor 2012a). We note that Stern & Laor found 
that EW(Ha) is independent of L at L < 10 46 erg s _1 , 
and therefore we assume that the same EW pertains at 
L 10 46 erg s _1 . The implied EW in velocity units is 
26,000 km s , and the continuum flux density is derived 
from this assumed EW. The units are chosen for easy 
conversion into the number of photons expected in an 
observation. 

Using eqn. (13) in eqn. (12), the implied S v for the 
simplified BLR is 


-F, 


lu<i;rot s i n * 




(jjcont 


TBLR COSJ- 


U rot Sin l 


, ( 14 ) 


which is shown in the bottom panel of Fig. 2, assuming 
j = 30° and 0(7 -blr) = 100 pas. The v = 5000 km s^ 1 
photons are offset from the v = —5000 km s _1 by « 
2r- B LRCOsj • $„/($„ + 4>“ nt ) = 150 pas. 


3.2. Local Line Broadening 

We now relax the assumption that the local line broad¬ 
ening at each location in the disk is negligible. Therefore, 

$*(r, p') oc S(r — r B LR)e ^ ( sm v -„ rot iisi) 

(15) 

In Figure 3 we plot and S v for different <j/(v ra t sinz). 
In order to match the emission line-profile of a putative 
observation, u ro t is chosen so that the line width remains 
constant in the different models. For convenience, we 
characterize the line width by the Half Width at Half 
Maximum (HWHM). 

Figure 3 illustrates that as cr/(v I0 t sinz) increases, 
loses its double-peaked shape and the maximum spec- 
troastrometric signal IS^I decreases. This behavior oc¬ 
curs because as a/(v ro t sinz) increases, photons from a 
wider range of ip' contribute to a given value of v, thus 
smoothing the spectral profile and S v . 

What is the value of a expected in the BLR? There 
are some indirect observations which suggest that a ~ 

5 The value of 3>*(r, </?') is equal to the integrand in eqn. 7 in 
CH89, divided by the observed photon energy. 



Figure 3. The emission line profile (upper panel) and photocen¬ 
ter offsets (lower panel) of a BLR ring with local line broadening. 
At each location in the ring, the emitted photon wavelength distri¬ 
bution is assumed to be a Gaussian with width a. The black lines 
assume negligible broadening as in Fig. 2, while other plotted lines 
assume the noted value of cr/{v ro t sin i). As ry/(v T ot sin*) increases, 
the two peaks in the line profile merge into a single peak, and the 
maximum value of |Su| decreases. 

bVot sin i ~ 0.5u ro t 6 . In most quasars the BLR pro¬ 
files do not exhibit a double-peaked profile (Strateva 
et al. 2003), which is expected when cr < 0.5u ro tsint 
(see top panel of Fig. 3). This lack of a double-peak is 
sometimes attributed to the BLR being emitted from 
a wide range of r, however this explanation is disfa¬ 
vored by the RM results which suggest that the BLR 
originates from a narrow range in r (see §1). There¬ 
fore, the single-peaked profile typically observed is more 
likely due to er /v TO t sini of order unity. Furthermore, Os- 
terbrock (1978) argued that the statistics of BLR line 
widths suggest a/v TOt ~ 0.4, which is also roughly con¬ 
sistent with the inferred BLR covering factor (CF) of 
~ 0.3 (Korista et al. 1997), assuming that a/v TOt is 
a measure of the height-to-diameter ratio of the BLR 
disk. Due to the above arguments, for the purpose of 
estimating the strength of the spectroastrometric signal 
we will henceforth assume a/(v TO t sini) « 1. The bot¬ 
tom panel of Fig. 3 shows that <r/(n ro t sin?) « 1 implies 
SiHWHM ~ ±0.25 0(r B LR)- 

3.3. Distribution in r 

We now relax the idealized assumption that the BLR 
originates from a single r, and consider a general spatial 
distribution of the broad line emitting gas, parameterized 

6 The typical value of sin i is expected to be ^ 0.5 in the standard 
AGN unification model, where the broad line region is visibile only 
if the quasar is viewed at pole-on inclinations. 
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by /(r), the line emission per unit logr: 


$ 


*/ n f(r) -feotljn0i/ siB 

„(r,^) =-e ^ V 


rr ) 2 


( 16 ) 


What is the expected shape of f(r)7 As noted in §1, RM 
studies suggest that the BLR emission originates from a 
small dynamical range in r, i.e. f(r) is strongly peaked 
around r = trm- We therefore, parameterize f(r) as a 
double power law f(r) oc r ±a , where the index is posi¬ 
tive at r < trm and negative at r > trm, and calculate 
S v for a = 1 and 2. These f(r) are shown in the top 
panel of Figure 4, together with a (5-function (thin ring) 
for comparison with Figs. 2 and 3. We assume that /(r) 
covers the range 0.03 t-rm — 30trm- Additionally, we also 
consider a physically motivated /(r), where f(r) drops 
at r > 2 crm = r su b due to absorption of the ionizing 
photons by dust grains, and also drops at r < trm due 
to collisional suppression of the line emission. This f(r) 
was calculated for different lines by Baskin et al. (2014), 
assuming a constant BLR CF per unit logr. Since below 
we estimate S v for Ha, we use the Baskin et al. calcula¬ 
tion of f(r) of Ha, which is similar in shape to f(r) of 
H/3 seen in their fig. 5 7 . We refer the reader to Baskin et 
al. for the details of the calculation. 

The middle and lower panels of Fig. 4 show the im¬ 
plied and S v , for the different /(r) shown in the 
top panel. As above, we assume 0 (crm) = 100 /xas, 
i = j = 30°, an EW of a typical Ha (= 26 000 km s -1 ), 
and v rot at r = Trm such that HWHM = 5000 km s -1 . 
We assume v ro t (r) oc r -1 / 2 , and based on the argu¬ 
ments in the previous section, cr(r)/v TO t(r) sin i = 1 . 
The middle panel shows that 4>„ of the different f(r) 
are similar, up to somewhat more extended line wings 
in the /(r) oc r^ 1 and Baskin et al. distributions com¬ 
pared to the f(r) oc r ±2 and (5-function distributions. 
This similarity is expected since all assumed f(r) drop 
sharply when r is significantly larger than or significantly 
smaller than r^M■ In contrast, the shape of S v seen 
in the lower panel differ significantly for the different 
f(r). At v = 2000 km s -1 , for example, S v = 25/xas 
for f(r) oc r ±2 , compared to S v = 60/ms for f(r) oc r ^ 1 
and S v = 90/xas for the f(r) calculated by Baskin et 
al. However, at |u| ss HWHM = 5000 km s _1 , the 
dependence of S v on f(r) is significantly weaker than 
at v = 2000 km s^ 1 , where S„ = hwhm is in the range 
25 — 37 /ras for all considered f(r). 

Why is S v at v < HWHM so sensitive to emission from 
large r? As seen in eqn. ( 12 ), S v is a photon-weighted 
average of the contribution to S v from each r. Eqn. (14) 
demonstrates that the contribution to S v from a specific 
r scales as oc rv~^ t oc r 3 ^ 2 , for v < v TO t (r) sin i. While 
eqn. (14) is derived assuming er <C u ro tsini, the bottom 
panel of Fig. 3 shows that it is also roughly a correct 
description of S v for general er, up to multiplication by 
a constant which depends on er. Therefore, in the case 
where /(r) oc r _1 at r > trm, this r 3 / 2 term implies that 
S v is dominated by photons from the largest r with sig¬ 
nificant emission at v, despite that these large-r photons 
have a negligible contribution to the line profile. This 


7 Baskin et al. plot f(r) in terms of the line EW, assuming a CF 
per unit log r of 0.3. We use their BLR model with solar metallicity 
and an ionizing spectral slope of —1.6. 
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Figure 4. (Top) Possible radial distributions of the line emission. 
The black arrow is a 5-function (thin ring), as used in Figs. 2 and 3. 
The red solid line and the green dashed line are parameterizations 
of f(r) which are consistent with RM studies, where f(r) peaks at 
r = r^M and drops at higher and lower r. The top axis notes the 
physical scale for a quasar with T 1450j ^ = 10 47 erg s -1 . The blue 
dash-dotted line is a physically-motivated Ho-emission distribution 
from Baskin et al. (2014, fig. 5 there), where the line emission drops 
at r > r su b = 2r^M due to absorption of the incident ionizing pho¬ 
tons by dust grains, and also drops at r <C due to collisional 
suppression of Ha. (Middle) The line spectral profiles for the dif¬ 
ferent f(r) shown in the top panel, assuming cr/(v TO t sini) = 1 and 
HWHM = 5000 km s —1 . The spectral profiles are similar to each 
other. (Bottom) The expected photocenter offsets for the differ¬ 
ent f(r). At |u| < HWHM, S v is highly sensitive to the relatively 
weak emission from r trm- At |u| « HWHM, S v only weakly 
depends on the assumed f(r). 

effect is even more pronounced for the physical f{r) cal¬ 
culated by Baskin et al. (2014), in which more photons 
originate from r > 10 trm than in the f(r) oc r^ 1 case. 

To conclude, the spectroastrometry signal at v < 
HWHM is very sensitive to f(r) at r > trm, and there¬ 
fore can be used to probe the line emission from this 
regime, as further discussed in § 6 . At v comparable to 
the HWHM, S v probes trm, be. S v probes the scale at 
which the bulk of the BLR emission is produced. 

3.4. Contaminants 
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In a real observation, some of the photons with wave¬ 
lengths spanned by the observed broad line originate 
from non-BLR sources, and therefore we need to account 
for their affect on S v . We have already referred to the 
effect of quasar continuum photons, which dilute S v . and 
can be accounted for by dividing the observed S v by the 
fraction of line photons at each v (see §3). In this sec¬ 
tion, we discuss additional non-BLR contaminants and 
how to account for their effect on S v . 

3.4.1. Narrow emission lines 

Photons from the NLR originate from r tblr, and 
hence any asymmetries on NLR scales (as seen in fig. 5 
of Bailey 1998) may dominate S v at v with narrow line 
emission. The largest relevant r is rpsF, the physical 
size spanned by the PSF, since emission from r > 7 "psf 
can likely be filtered out during the data reduction phase 
of the spectroastrometric observation. For an assumed 
PSF = 70 mas we get tpsf « 10 3 tblr (^psf = 600 pc at 
z = 2). Now, since in the calculation of S v (eqn. 12) there 
is a factor of r in the integrand, NLR emission which 
originates from IO'Wblr may dominate S Vl even if the 
NLR flux density is only ~ 0.1% of the broad line flux 
density. This is the same argument made in §3.3 for why 
weak BLR emission from large scales can dominate S v . 
Hence, photons from wavelengths near narrow emission 
lines will need to be disregarded when measuring the 
photocenter of the BLR. 

For example, at a broad Ha luminosity of 10 44 erg s _1 , 
which is equivalent to L = 10 46 erg s _1 (Stern & Laor 
2012a), the mean fluxes of the narrow Ha, [O i], [S il] 
and [N ii] lines are 0.3 — 3% of the broad Ha flux (fig. 3 
in Stern & Laor 2012b, fig. 1 in Stern & Laor 2013). The 
flux density ratio of [N ii] and the narrow Ha is ~ 10 
times this value for a typical NLR-to-BLR FWHM ratio 
of 0.1, while [O i] and [S Ii], which reside on the broad Ha 
wings, have even higher flux density ratios. Therefore, 
since the NLR to BLR flux density ratios is 0.1% at v 
with these strong narrow lines, the relevant S v will likely 
be dominated by the NLR. 

We search also for weaker lines in the CLOUDY (Ferland 
et al. 2013) model of the NLR described in Stern et al. 
(2014) 8 . We find three additional lines which may be 
strong enough to dominate S v . The lines are [S ill], Ar V, 
and He I, at = —11 500, 5840, and 5270 km s _1 , and 
have fluxes which are 0.08, 0.03, and 0.01 times the flux 
of the narrow Ha, respectively. The latter two lines can 
be seen as small bumps in the high resolution spectra 
of NGC 4151 (fig. 1 in Arav et al. 1998). Photons near 
these lines may also have to be disregarded. In §5 below 
we demonstrate that the emission from even weaker lines 
will not significantly affect S v . 

We note that the above estimates for the NLR strength 
is likely an upper limit, since spectroastrometry candi¬ 
dates have L 10 46 erg s _1 , while the narrow to broad 
flux ratio decreases with L as ~ L~ 03 , probably due to 
the decrease in NLR covering factor with increasing L 
(Stern & Laor 2012b). Additionally, some of the narrow 
line emission may originate from r > rpsF and therefore 

8 We use the Stern et al. (2014) NLR model with an ionizing 

1 /2 

slope of — 1 . 6 , solar metallicity, and r = 600 pc, which corre¬ 
sponds to r — 1000 tblr ~ ^PSF 5 where the NLR contribution to 
S v peaks. 


can be filtered out. 

3.4.2. Continuum emission from the host galaxy, distant 
scattering media, and dust grains 

At optical wavelengths, the polarized flux density in 
unobscured AGN is typically 0.5 — 5% of the total flux 
density (Smith et al. 2002), suggesting that some of the 
observed emission is scattered emission. The scattering 
medium may reside on scales tblr, and therefore, 
based on the same reasoning used above for the NLR, 
asymmetries in this scattering medium can in princi¬ 
ple dominate S v . However, because a distant scatter 
views the BLR as a point source, the opening angle 
covered by the scatterer will not be a function of lo¬ 
cation in the BLR, and hence not a function of veloc¬ 
ity. Since the scattering efficiency is a weak function of 
A, and a v = 10 000 km s _1 emission line spans merely 
dA/A = 0.03, the scattering efficiency is also a weak func¬ 
tion of velocity. Therefore, we expect the contribution of 
the scattered flux to S v to be roughly constant across 
the emission line and the flanking continuum. We can 
account for the effect of a distant scatterer by interpolat¬ 
ing S v from the continuum flanking the broad line profile, 
and subtracting out the result. This approach will also 
account for the effect of the host galaxy emission on S v , 
since any asymmetry in the host galaxy is also likely to 
be a weak function A. 

We note that the rotation of the polarization angle 
across the Ha profile found by Smith et al. (2002, 2005) 
suggests the existence of scattering media also at r > 
i’blr- Such nearby scatterers will not dominate S v . 

If one measures the spectroastrometric signal of a rest- 
frame IR emission line such as Paa, then the dust IR 
continuum can affect S v . The S v profile of the dust con¬ 
tinuum is also expected to be a weak function of A, and 
therefore can be accounted for using the same interpo¬ 
lation method noted above for accounting for scattering 
media. However, since the dust continuum is a signifi¬ 
cant fraction of the flux at IR wavelengths, in contrast 
with the host and scattered flux, any second-order de¬ 
pendence of S v on A which is not accounted for may 
significantly affect S v , and therefore may pose a problem 
for performing spectroastrometric on rest-frame IR lines. 

4. OBSERVATIONAL CONSIDERATIONS 
4.1. Systematics 

We note that most of the atmospheric and instrumen¬ 
tal systematic errors relevant to imaging astrometry are 
irrelevant to BLR spectroastrometry, since in BLR spec¬ 
troastrometry one measures the astrometric position of 
the broad emission line relative to the adjacent contin¬ 
uum, rather than in any absolute sense. Specifically, 
the red and blue sides of the broad emission line are 
expected to be offset from the continuum in different di¬ 
rections (Figs. 2-4). Most systematic effects on the 
photon centroid are expected to be a weak function of A, 
since the line and continuum photons travel in the ex¬ 
act same way through the atmosphere and instrument, 
and are only split in the wavelength domain. Therefore, 
to mimic a spectroastrometric signal, any instrumental 
systematic needs to vary so strongly with A that within 
dA/A « 0.03 it produces a relative offset between the 
red and blue wings of the emission line. Any systematic 


which is a weak function of A can be accounted for by 
interpolating the centroid of the continuum flanking the 
emission line, a method which accounts also for the ef¬ 
fect of astrophysical continuum sources such as the host 
galaxy or scattering media (§3.4.2). The error of this 
interpolation is limited by the photon count of the con¬ 
tinuum, and is estimated below. 

On the other hand, instrumental systematics which are 
a strong function of A may be a limiting factor, rather 
than the number of photons collected. Pontoppidan et al. 
( 2011 ) showed that rotating the slit by 180 deg can mit¬ 
igate such systematics which are related to flat-fielding 
and removal of telluric lines. This method is based on 
the fact that the spectroastrometric offset is a real offset 
between line photons on the sky, and therefore a 180 deg 
rotation of the slit flips the direction of the spectroas- 
trometry offset in detector coordinates, while most in¬ 
strumental systematics would be fixed in detector coor¬ 
dinates. By applying this 180 deg flip, Pontoppidan et 
al. reached photon-limited accuracies of 100 — 500 /xas, 
i.e. at an accuracy of 100 pas their observations were still 
free of systematic effects. We therefore find it reasonable 
that systematics are, or can be reduced to ~ 25 pas, the 
accuracy level required for BLR spectroastrometry. 

An additional method to differentiate between a true 
signal and an instrumental systematic would be to test 
the predicted j-dependence of the spectroastrometric sig¬ 
nal (eqn. 12), using multiple slit orientations. In order 
for an instrumental systematic to mimic the spectroas¬ 
trometric signal it would to have to behave exactly the 
same way with rotation of the slit. 

4.2. Choice of Target and Emission Line 

The choice of an optimal target and emission line for 
spectroastrometry is motivated by four main factors. 
The angular size of the emitting region, the photon flux of 
the line, the strength of the non-BLR contaminants listed 
in the previous section, and the observed wavelength of 
the line. The Ha line scores highly given all of these 
considerations, since it is the strongest optical emission 
line, it is emitted from the outer part of the BLR (Bentz 
et al. 2010 ), and its wavelength is near 1 pm where the 
accretion disk and dust continuum emission are at a min¬ 
imum. Furthermore, at z = 2 — 2.5 where 0(rRM) is near 
its peak value (Fig. 1), Ha is redshifted into the near-IR 
K-band where adaptive-optics works best. 

Also, we note that the idea of an ordered rotational 
velocity field in the BLR is best established for low- 
ionization lines, such as the Balmer lines. Higher ioniza¬ 
tion lines appear also to have an outflowing wind com¬ 
ponent (Richards et al. 2011, see further discussion in 
§ 6 ). Hence, interpreting the spectroastrometic signal of 
high-ionization lines will likely be more complicated than 
interpreting the signal of low ionization lines. 

Given the above arguments, below we simulate the ex¬ 
pected spectroastrometry signal for Ha. We consider also 
the prospect of measuring S v on Mg II, which falls in the 
NIR bands at z = 3 — 7.5, and also S v of Paa in the local 
quasar 3C 273, the brightest quasar on the sky. 

4.3. Choosing the Slit Orientation 

Since S v oc cos j (eqn. 12 ), the required number of 
photons to achieve a given S/N scales as cos ” 2 j, where 


j is not known a-priori. In the context of measur¬ 
ing the spectroastrometric signal of stellar disks, Pon¬ 
toppidan et al. (2008) suggested an observing strat¬ 
egy where one observes the target at three slit ori¬ 
entations rotated by 60° from one another. Since 
(cos 2 (j) + cos 2 (j — 60) + cos 2 (j — 120)) /3 = 0.5 for all 
j, than the required observing time with this technique 
is twice the observing time required for the j ' = 0 case. 
Therefore, below we assume j = 0, and multiply the re¬ 
quired observing time by a factor of two to account for 
the unknown j. 

We note that even if one performs a spectroastromet¬ 
ric observation with an Integral Field Unit (IFU), which 
maps the sky into 2D pixels and preforms simulatenous 
spectroscopy in each pixel, rotation of the unit is still 
likely to be required. An IFU maps one spatial direc¬ 
tion continuously onto the detector, while the other spa¬ 
tial direction is mapped discontinuously. Constructing 
the spectroastrometic signal along the spatial dimension 
which is mapped continuously onto the detector is tan¬ 
tamount to constructing the signal in a long-slit obser¬ 
vation. However, constructing the spectroastrometic sig¬ 
nal along the dimension which is mapped discontinously 
requires modelling the projection of this spatial dimen¬ 
sion onto the detector to the level of 10” 3 of a pixel, 
as required for BLR spectroastrometry, and is therefore 
unlikely to be feasible. Therefore, using an IFU has no 
clear advantage over using a long-slit. 

5. SIMULATED SPECTROASTROMETRY SIGNAL 

In this section we simulate the spectroastrometric sig¬ 
nal for a z = 2 target on 8 m and 39m telescopes, and then 
consider measurements of targets at other redshifts. 

5.1. The Ha Line of SDSS J1521+5202 

SDSS J152156.48+520238.5 (hereafter J1521+5202) is 
the most luminous SDSS quasar where Ha falls in the Ri¬ 
band (Fig. 1). J1521+5202 has L 1450 a = 10 476 erg s” 1 , 
0(rRM) = 111/zas, and 2 = 2.208. The FWHM of its 
Mg II line is 9300 km s” 1 (Shen et al. 2011). As the NIR 
spectrum of J1521+5202 is not available, we use as a tem¬ 
plate the spectrum of SDSS J154019.57-020505.4 (here¬ 
after J1540-0205), a L 1450 a = 10 45 ’ 4 erg s” 1 quasar at 
2 = 0.3 with FWHM(Ha) = 9730 km s” 1 . The two spec¬ 
tra should be similar since the FWHM of the Balmer lines 
is correlated with the FWHM of Mg II (Shen et al. 2008). 
We multiply the observed photon flux of J1540-0205 by 
the photon flux ratio of the two quasars, 10 47 ' 6 /10 45 ' 4 ■ 
(d L (z = 2.2 )/d L (z = 0.3)) 2 • (1 + 2.2)/(l + 0.3) = 3.1. 
The normalized spectrum is shown in the top panel of 
Figure 5, together with the broad Ha, narrow lines and 
continuum fits described in Stern & Laor (2012a), which 
we utilize below. We note that since J1521+5202 is more 
luminous than J1540-0205, the host galaxy is expected 
to be relatively weaker, as are the narrow emission lines 
(Stern & Laor 2012b). We determine the flux of the weak 
[S in], Ar V, and He I lines from the flux of the narrow 
Ha and the model flux ratios mentioned in §3.4.1. These 
lines are too weak to be observable in the SDSS spec¬ 
trum. 

In order to estimate S v , the broad emission line is 
modeled as a rotating disk with cr/(u ro t sin*) = 1 and 
either /(r) oc r ±2 , or the f(r) based on Baskin et 
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Figure 5. The simulated spectroastrometric signal of the broad 
Ha in J1521+5202, a T 1450 ^ = 10 47,6 erg s —1 quasar at z = 2.2. 
(Top panel) The yellow line plots the normalized spectrum of 
J1540-0205, used as a template for J1521+5202. The black line 
plots the fit to the continuum, broad Ha, and the noted narrow 
lines (marked by vertical dotted lines). (2 nd panel) The expected 
S v , assuming cr/(u ro t sini) = 1 and the f(r) oc r± 2 shown in Fig. 4. 
The narrow lines likely have asymmetry on large scales, therefore 
at v with strong NLR emission S v is on the scale of the PSF, which 
is beyond the plotted axis. The weak Ar v and He I lines have a 
comparable contribution to S v as the broad Ha. At v which are not 
dominated by narrow lines (thick line), the broad Ha and quasar 
continuum dominate S v . The photocenter at v = 4000 km s —1 
is offset from the photocenter at v = —4000 km s —1 by 85/xas. 
(3 rd panel) The solid line marks the expected S v at v where the 
NLR does not dominate S v . The error bars are the simulated 
centroid statistical uncertainties, for spectral bin sizes of 1000 and 
4000 km s —1 in the emission line and continuum, respectively. We 
assume a 10 hr integration on an AO-assisted 8m telescope, with 
Strehl = 0.4 and /i ns t = 0.2. (Bottom panel) The expected S v 
and centroid uncertainties assuming the physically-motivated f(r) 
shown in Fig. 4. 


al. (2014). We assume for simplicity that the angu¬ 
lar size of the Ha-emitting region is #(rRM), though 
the true 9 may be higher since trm was determined on 
H/3, and Ha is sometimes observed to have longer RM 
lags than H/? (Bentz et al. 2010). We set i> ro t(ntM) sini 
of the model such that the FWHM in the model pro¬ 
file equals the observed FWHM of 9730 km s^ 1 , i.e. 
2u ro t (t’rm ) sin i = 7400 km s -1 for the /(r) oc r ±2 case 
and 2i! rot (rRM) sini = 11000 km s -1 for the f(r) from 
Baskin et al. We normalize the model <E» (= / ^du) to 
give the same value of the observed $. We note that at 
|v| < 5000 km s _1 , the observed profile of the broad 
Ha differs by up to 20% from the <!>„ calculated by the 
BLR model. However, since the goal of this study is 
to estimate S v rather than provide a model which accu¬ 
rately fits the observed emission line, we disregard this 
difference. 

In order to estimate S v of the narrow lines, we as¬ 
sume the worst-case scenario where the narrow lines orig¬ 
inate from r = rpsF = 600 pc (assuming an 8m-telescope 
diffraction limited PSF of 70mas), where the contribu¬ 
tion to S v is maximal (§3.4.1). For simplicity, we assume 
the NLR emission also originates from a turbulent rotat¬ 
ing disk with r> rot sin i chosen to fit the observed narrow 
lines FWHM of 450 km s -1 . The total photon flux of the 
NLR model is chosen to fit the observed photon fluxes 
of the strong emission lines (Ha, [O i] [N II], and [S n]), 
and the model photon fluxes of the weak emission lines 
([S III], Ar V, He II). Finally, quasar continuum photons 
are assumed to originate from r = 0. 

The second panel in Fig. 4 shows S v versus v for f(r) oc 
r ±2 . Thin lines mark v values which are contaminated by 
emission from narrow lines. The NLR spectroastrometric 
signal is on the scale of the PSF, which is beyond the 
plotted axis. Note that the weak signal of He I justifies 
our choice in §3.4.1 to disregard narrow lines with a flux 
which is lower than the flux of He I. The thick line marks 
the range of v where the broad Ha and quasar continuum 
dominate S v , and only these v are used further in the 
analysis. The photocenters at the blue and red wings 
are offset by up to 85 /zas from each other. 

To estimate the expected error on the measured photo¬ 
center, we assume a t = 10 hr observation on an 8m tele¬ 
scope. Using the specifications of the Nasmyth Adaptive 
Optics System (NAOS) and COude Near Infrared CAm- 
era (CONICA) as guidance 9 , we assume a collecting area 
of A = 38m 2 , a diffraction-limited FWHMpsf = 70 mas, 
a Strehl ratio of 0.4, and an instrument photon collect¬ 
ing efficiency of /; nst = 0.2. For the observed ~ 

10 6 [m _2 hr _1 (l0 3 km s -1 ) 1 ] at half maximum of the 
broad Ha (top panel of Fig. 4), the implied number of 
collected photons per 1000 km s -1 is 

—= $> v -A-t-Strehl-fi ns f0-5 = 15xl0 6 (l0 3 kms _1 ) 1 , 

du 

. T) 

where the factor of 0.5 reduction arises from the slit ori¬ 
entation as discussed in §4.3. 

The value of dA r ph/d'c computed in eqn. (17) implies a 


9 See user manual at http://www.eso.org . 
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Figure 6. Simulated spectroastrometric signal with a 39m tele¬ 
scope. As in the bottom panels of Fig. 5, the solid line denotes 
the expected S v at velocities in which the BLR and quasar con¬ 
tinuum dominate S v • Error bars denote the expected statistical 
uncertainty in the photocentroid measurement. In both panels, 
f(r) from Baskin et al. (2014) is assumed. (Top panel) The sig¬ 
nal expected on the broad Ha of J1521+5202, the object which 
appears also in Fig. 5. The higher photon flux and smaller PSF in 
39m telescopes compared to 8m telescopes enables constraining S v 
to an order of magnitude higher velocity resolution, with a frac¬ 
tion of the observing time. (Bottom panel) The signal expected 
on the broad Mg II line of SDSS J114816.64+525150.3, a, z = 6.4 
quasar. 


v < HWHM the maximum expected S v is a factor of 
four larger than when assuming f(r) oc r ±2 , and the rel¬ 
ative errorbar sizes are correspondingly smaller. 

In the top panel of Figure 6, we show the expected 
spectroastrometric signal of J1521+5202 for a 10 min ob¬ 
servation on a next-generation 39 m telescope, zoomed-in 
on the emission line for clarity. We assume that A oc d 2 , 
FWHMpsf cx d _1 , a Strehloi 0.4 (e.g. Clenet et al. 2010), 
and /i ns t = 0.4 (e.g. Davies et al. 2010). Hence, eqs. 17 
and 18 imply that e oc d~ 2 , and that the required observ¬ 
ing time to reach a certain value of e scales as d~ 4 . We 
therefore decrease the spectral bin sizes to 200 km s -1 . 
With a next generation telescope, one can derive S v to a 
high velocity resolution, and therefore yield strong con¬ 
straints on the BLR kinematics. 

5.2. Spectroastrometry Estimates for Additional 
Quasars 

In this section, we estimate the statistical S/N of the 
spectroastrometric signal of the most luminous quasars 
at different 2. To this end, we define S/N(S' re d — SWe), 
to be the photocenter offset between the red and blue 
wings of a broad emission line, which is an aggregate 
S/N that effectively combines the data points shown in 
Figs. 5 and 6. For each wing, we include all photons 
with velocities 0 < |u| < FWHM, excluding velocities 
which might be contaminated by narrow line emission. 
We address quasars in which Pact, Ha, or Mg II fall in 
one of the NIR transmission windows, and require A res t 
is at least 10 000 km s -1 away from the edge of the at¬ 
mospheric window, in order to adequately sample the 
quasar continuum and avoid significant telluric absorp¬ 
tion. These quasars are listed in Table 1, which appears 
in the Appendix. The S/N of less luminous quasars can 
then be estimated by noting that 

+rm) L 1/2 

b/N(z) oc-oc- —pr oc L (19) 

N~ 1/2 


statistical la photocenter positioning error e v of 

FWHMpsf 


2.35 (^) 


1/2 


= 7.7 pas 


(18) 


The third panel in Fig. 5 shows the expected errors in 
the photocenter measurements, for bins of size Av = 
1000 km s -1 in the red and blue wings, and bins of size 
Av = 4000 km s^ 1 for the red and blue continuum re¬ 
gions. To simulate an observation, the center of the sim¬ 
ulated errorbars is offset from the expected signal (solid 
line) by a random offset chosen from a Gaussian dis¬ 
tribution with dispersion e v . Fig. 5 demonstrates that if 
systematics can be reduced to the level of a few times the 
value of e v calculated in eq. 18, then the photocenter off¬ 
set of the red wing from the blue wing is detectable with 
existing telescopes. For comparison, Pontoppidan et al. 
(2011) reached photon-limited uncertainties of 100 pas in 
the context of stellar disks, but this precision level was 
set by the number of photons collected rather than by 
systematic errors. 

In the bottom panel of Fig. 5, we show the expected 
spectroastrometric signal assuming the /(r) based on 
Baskin et al. (2014). As expected from Fig. 4, at 


where we used eqs. 1 and 18 for the dependence of 0(?~rm) 
and e v on L, respectively, and emphasize that eqn. (19) 
is only valid at a fixed redshift. We assume a Strehl 
ratio of 0.2, 0.4, and 0.4, for the J, H and K bands, 
respectively, and /i nst = 0.2 (0.4) for all bands for an 8m 
(39m) telescope 10 . 

Potential narrow line contamination near Mg n and 
Paa is derived in a similar fashion as for narrow lines near 
Ha (§3.4.1). Narrow lines with a large enough flux to sig¬ 
nificantly effect S v near Mg II include the narrow Mg II 
doublet, a He II line at UMg n = —7000 km s _1 , an Ar IV 
line at VMg u = 6000 km s+ and a Mg V line at «M g n = 
6100 km s _1 . Similarly, near Paa there are the Br<5 and 
Bre lines, three He lines at vp aa = —255, —1055 and 
—1855 km s _1 , and an [Fe n] line at vp aa = 2140 km s _1 . 
All photons with v within 750 km s _1 of these lines are 
excluded from the analysis. We note that some of the 
mentioned lines are permitted lines, and therefore should 
also exhibit broad line emission, however these broad 
lines are likely too weak to affect S v significantly. 

10 Instrument throughput based on NACO and MICADO 
(Multi-AO Imaging Camera for Deep Observations, Davies et al. 
2010) specifications. 
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The rotating disk parameters of the model are the same 
as in the previous section, with u ro t(rRM) and the nor¬ 
malization chosen so the model reproduces the estimated 
values of FWHM and $ listed in Table 1 in the Appendix. 
Since not all targets have available NIR spectra, we es¬ 
timate FWHM and $ from observed emission lines, as 
detailed in the Appendix. For Mg II, we assume that 
the total flux is contributed from the two lines in the 
doublet (separated by 770 km s _1 ) according to a 2 : 1 
ratio. Therefore, at a given v the expected centroid is 
the photon-flux-weighted centroid of the individual lines 
in the doublet. 

The value of S/N(5 , re d — Sbiue) is the number of stan¬ 
dard deviations from which Sped — /Sfoue deviates from a 
straight line interpolation of the photon centroids of the 
continuum bins. It is equal to 

S/N(-SVed - Sblue) = fl-ed-Sb! (2Q) 

( £ red + £ blue + £ cont) 

where £bi U e and e re d are the statistical errors on the pho¬ 
tocenters of blue and red wing photons, respectively, cal¬ 
culated from the expressions in eqs. 17 and 18. The error 
on the continuum level, e C ont, is equal to the error on the 
evaluation of a straight line interpolation from the data 
points in the blue and red continuum regions: 

I^red ^blue [2 i 2 

Ccont “ _ ( £ red cont ' £ blue cont ) > 

^red cont <T>lue cont 

(21) 

where Cbiuecont and e re dcont are the statistical errors on 
the photocenters of each continuum bin (FWHM < |u| < 
2 FWHM), and the continuum flux is calculated from the 
emission line EW noted in Table 1. The values of n re d and 
fbiue are the photon-flux average v of the red and blue 
wing bins, respectively, while i; r edcont = -i+iuecont = 
1.5 FWHM are the average velocities of the continuum 
bins. We note that in the S/N calculations shown in 

Fig- 7 , £ cont/( £ red + e biue) = °- 07 - °- 3 ; and therefore 
including e CO nt in the calculation does not significantly 
affect S/N^SVed — ■S'biue)- However, in a real observa¬ 
tion the amount of continuum photons may be limited 
if the emission line falls near the edge of an atmospheric 
window, which will increase e con t, and hence decrease 
S/N(S red - Sbhie). 

The expected S/N(5' re d — S'biue), for 8m and 39m tele¬ 
scopes, are shown in Figure 7. At 1 < z < 2.5, a night 
on an 8m is sufficient to achieve S/N(S re d — <5biue) ~ 10, 
even when assuming /(r) oc r~ 2 . Assuming the Baskin 
et al.-based /(r) increases the expected S/N by a factor 
of w 2, while an hour on a 39m telescope increases the 
expected S/N by a factor of ~ 10 compared to 10 hours 
on an 8m. Fig. 7 suggests that with next-generation tele¬ 
scopes, the spectroastrometic signal can be detected for 
quasars as distant as z = 6.5. As an example, in the bot¬ 
tom panel of Fig. 6 we plot a simulated observation of 
the Mg II line in SDSS ,1114816.64+525150.3, a z = 6.43 
quasar discovered by Fan et al. (2003). 

6 . DISCUSSION 

6.1. The Assumption that the BLR is a Rotating Disk 

For the purpose of considering ordered motion, we as¬ 
sumed that the BLR has a flattened geometry and that 
its kinematics are dominated by rotation. These assump- 
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Figure 7. Expected S/N for spectroastrometry candidates at dif¬ 
ferent z. Candidates are chosen from the most luminous quasars at 
each 2 : (Fig. 1), where a strong low-ionization emission line falls in 
one of the NIR bands. The y-axis is the S/N of the centroid offset 
of the blue and red wings of the emission line (eqn. 20), for 10 hr 
on an 8m telescope (small symbols), and 1 hr on a 39m telescope 
(large symbols). The marker type denotes the assumed radial dis¬ 
tribution of the line emission (Fig. 4), and bold markers denote the 
examples shown in Figs. 5-6. The spectroastrometric signal of 
the BLR in 1 < z < 2.5 quasars is detectable at S/N > 10 with 
existing telescopes, while next-generation telescopes will be able to 
observe the spectroastrometric signal up to 2 : = 6.5. 
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tions are supported by several observations of the low- 
ionization BLR gas, which we briefly summarize below. 

Evidence for a flattened geometry comes from the AGN 
unification picture. According to the the standard AGN 
paradigm, in unobscured AGN (known as type 1 AGN, 
including quasars) the accretion disk is viewed at an¬ 
gles near face-on. When an AGN disk is viewed edge-on, 
the accretion disk and broad line emission are obscured 
by a geometrically thick dusty ‘torus’ (Antonucci 1993). 
Maiolino et al. (2001) and Gaskcll (2009) argued that 
if the BLR is spherically symmetric, then a fraction of 
the quasars, equal to the BLR covering factor of ~ 30%, 
should be viewed through the BLR clouds. Such BLR 
gas along the sightline would produce Lyman-edge ab¬ 
sorption and low ionization broad absorption lines (LoB- 
ALs). However, only 1% of quasars are LoBALs (Trump 
et al. 2006), in contrast with the prediction of a spheri¬ 
cally symmetric BLR. Hence, a more likely picture is that 
the low-ionization BLR resides predominantly near the 
disk-plane in a flattened geometry, along line of sights 
where the quasar would be viewed as obscured. 

As for our assumption of a rotating BLR, rotation is 
preferred over an outflow-dominated BLR, since in a sce¬ 
nario where the BLR is an outflow off the face of the 
disk, the redshifted and blueshifted photons will origi¬ 
nate from opposite locations along the minor-axis of the 
projected BLR ellipse, in contrast to along the major- 
axis in the rotating BLR scenario (see top panel of Fig. 
2). This geometry implies a shorter response time of the 
blue wing compared to the response time of the red wing, 
which is ruled out in velocity-mapped RM studies of low 
ionization lines (Maoz et al. 1991; Grier et al. 2013; Pan¬ 
coast et al. 2013), albeit for low-z low-L AGN. Further 
support for the rotating BLR assumption comes from 
the mean redshift of the broad H/3 relative to systemic 
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in SDSS quasars, which is consistent with the relativis¬ 
tic redshift expected in a rotating BLR (Tremaine et al. 
2014), and from the analysis of the microlensing signal 
of the quadruply lensed quasar HE 0435-1223 (Braibant 
et al. 2014). 

However, the high-ionization broad lines likely have 
a significant outflowing component. Evidence for an 
outflow can be seen in the line profile of C iv A1450 
(Baskin & Laor 2005; Richards et al. 2011), and in the 
high-ionization states typically observed in broad ab¬ 
sorption lines (Hamann 1997), which are predominantly 
blueshifted, suggesting an outflow. Since in an out¬ 
flow scenario the spatial separation of the redshifted and 
blueshifted gas is along the minor axis, which is smaller 
than the major axis by a factor of cosi, the expected 
S v should also be correspondingly smaller (assuming the 
same telr)- Since this signal reduction is relatively 
small, the spectroastrometric signal should be detectable 
also in an outflowing BLR. 

Note that as the rotational motion in the BLR de¬ 
creases relative to the random motion, S v is reduced as 
shown in the bottom panel of Fig. 3. In J1521+5202, for 
example, the S/N(S re d — Sbiue) = 17 expected for 10 hr on 
an 8m telescope (Fig. 7) is reduced to S/N(S re d — Sbiue) < 
5 if a /(u ro t sin*) > 3. Therefore, even a null detection of 
the spectroastrometric signal would place an interesting 
upper limit on the fraction of the BLR motion that is 
ordered. 

6.2 . Mbh Estimates 

Above we have demonstrated that by measuring the 
BLR photocenter offset S v , one can derive tblr- Given 
a measurement of tblr and an estimate of the Keple- 
rian velocity vk based on the line width, one can derive 
Mbh = rELR^K/G, analogous to RM-studies. Therefore, 
since any uncertainty in the tblr measurement propa¬ 
gates to the Mbh estimate, it is interesting to understand 
the implied uncertainty in 7 ~blr given a measurement of 

S v . 

For a given tblr, different values of a imply changes 
of a factor of two in S v (Fig. 3). Also, Fig. 4 shows that 
while different radial distributions of the line emission 
can change the S v profile significantly, the range of S v 
at |u| = HWHM for different f(r) spans merely a factor 
of 1.5. So, in the worst case scenario where nothing is 
known about a and /(r), for a given measurement of S v 
at | v | = HWHM one can estimate tblr to within a factor 
of ~ 3. However, this uncertainty can be reduced by con¬ 
straining a from the line profile, and by constraining f(r) 
using the measurements of S v at |u| < HWHM. For the 
purpose of estimating Mbh, these uncertainties on the 
t’blr estimate should be added to the uncertainty on vk 
due to the unknown inclination i of the disk plane to the 
line of sight. The total implied accuracy to which Mbh 
can be derived is an important topic for future study. 

Furthermore, er and f(r) are likely to be similar 
or the same for all BLRs. In this case, one can 
tie the spectroastrometry measurements of the BLR 
to lower-luminosity sources which have been reverbera¬ 
tion mapped, potentially removing this source of ‘noise’. 
With next-generation telescopes, one can potentially per¬ 
form spectroastrometry on objects which are low lumi¬ 
nosity enough to be also reverberation mapped. 


6.3. 8m-class Telescopes, 30m-class Telescopes and 

Space Telescopes 

Fig. 7 shows that existing telescopes can detect the 
BLR spectroastrometric signal in Ha in the most lumi¬ 
nous quasars at 1 < z < 2.5. Therefore, with existing 
telescopes one can use spectroastrometry to estimate the 
Mbh of the most luminous quasars at the peak of the 
quasar epoch. Also, the measured telr can be compared 
to the extrapolation of the tblr — L relation deduced by 
RM studies on lower luminosity AGN (eqn. 1). 

In §5 we show that the required integration time to 
reach a given S/N scales as d~ 4 , where d is the tele¬ 
scope diameter. Therefore, with next generation 30m- 
class telescopes, the required integration times are lower 
by a factor of (30/8) 4 = 200 than with 8m telescopes. 
Also, the smaller PSF of larger telescopes implies that 
reducing the systematics to a given angular precision is 
likely less challenging. 

The relatively short integration times required in 30m 
telescopes can be utilized to observe a large sample of 
quasars (from eqn. 20 and Fig. 7, reaching S/N = 10 
on quasars with L1450 = 10 47 erg s _1 at z « 2 requires 
7 — 15 minutes on a 39m). Thus, one can explore the de¬ 
pendence of S v on additional BLR parameters: e.g. the 
BLR FWHM, the EW of the lines, Mbh and L/LEdd, 
whether the high-ionization BLR shows signs of an out¬ 
flow (Richards et al. 2011), and whether the quasar spec¬ 
trum shows Broad Absorption Lines, which might be 
some proxy for orientation. The high photon flux can 
also be utilized to achieve a higher resolution in v, as 
shown in the top panel of Fig. 6, and thus further con¬ 
strain the kinematic model of the BLR. Additionally, 
with next generation telescopes one can explore fainter 
quasars, such as quasars as distant as z = 6.5 (Fig. 7), 
and fainter emission lines. This opens the possibility to 
perform spectroastrometry on multiple lines simultane¬ 
ously, at specific z where multiple lines fall in the NIR 
atmospheric windows (say C in] A1909 in J and Mg II in 
H at z ~ 5). Thus, one can constrain the BLR structure 
as a function of the ionization level of the line-emitting 

gas- 

Another interesting possibility is to perform BLR spec¬ 
troastrometry using the James Webb Space Telescope 
(JWST, Gardner et al. 2006), which has a diameter of 
6.6m and allows full coverage at 0.6 — 5 /rm without the 
issue of telluric absorption. Since Strehl = 1 in a space 
telescope, compared to Strehl = 0.4 assumed above for 
an 8m telescope, the observing time requirements are 
comparable ((6.5/8)~ 4 x 0.4 = 0.9). The continuous 
wavelength coverage can allow observing multiple emis¬ 
sion lines for the same quasar. Utilizing the Hubble 
Space Telescope (HST) for BLR spectroastrometry may 
also be possible, since the increased observing time due to 
the relatively small telescope diameter of 2.4m is partially 
offset by the lower diffraction limit in the UV. We differ 
a more thorough analysis of the possibility to perform 
BLR spectroastrometry with space telescopes to future 
work. 

6.4. Comparison with the Interferometry Results of 

Petrov et al. (2012) 

One can also use differential interferometry to measure 
the photocenter offset of the BLR relative to the contin- 
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uum. Petrov et al. (2012) applied this approach to the 
Pa line of 3C 273, and measured an angular radius of 
> 400 /zas. As noted there, this angular radius is larger 
by a factor of > 3 than the radius implied by RM of 
Ha. We suspect that this apparent discrepancy between 
fBLR measured by interferometry and tblr measured by 
RM is due to the high sensitivity of the photocenter to 
emission from large r (Fig. 4 and §3.3). Fig. 4 demon¬ 
strates that when a small fraction of the line photons 
come from r trm (he. when f(r) oc or when us¬ 
ing the physically-motivated f(r) calculated by Baskin 
et al. 2014), then the average S v at 0 < v < HWHM in¬ 
creases by a factor of 3 — 4 compared to when assuming a 
^-function radial distribution. This increase is consistent 
with the discrepancy found by Petrov et al. 

We note that a key advantage of spectroastrometry 
over interferometry is sky coverage. Currently interfer¬ 
ometers require co-phasing of the telescopes with a bright 
star or AGN (< 12 mag in the NIR), which needs to be at 
a small angular separation on the order of the isoplanatic 
angle at the science wavelength (typically < 20", Espos¬ 
ito et al. 2000). Only a handful of nearby quasars (Fig. 
1 ) are bright enough such that a nearby star is unnec¬ 
essary. In contrast, spectroastrometry can benefit from 
the ~ 50% sky coverage provided by Laser-Guide-Star 
Adaptive Optics (e.g. Diolaiti 2010). 

6.5. Constraints on Torus Models 

As first noted by Netzer & Laor (1993), the line emis¬ 
sion from gas at r > r su b which is exposed to the quasar 
ionizing radiation is suppressed, due to absorption of the 
ionizing photons by dust grains, rather than absorption 
of the ionizing photons by gas particles. This suppres¬ 
sion, seen as a drop at r > 2crm « r su b in the Baskin 
et al. 2014 /(r) profile shown in the top panel of Fig. 4, 
explains the apparent gap in line emission between the 
BLR and the NLR. This dusty gas at r > r su b, usu¬ 
ally referred to as the ‘torus’, also creates the dichotomy 
of obscured / unobscured AGN invoked in the standard 
AGN unification model discussed above. Due to the weak 
line emission of the dusty gas, the physical properties of 
the torus have hitherto mainly been constrained via the 
observed IR emission and via the demographics of ob¬ 
scured and unobscured AGN. However, in Figs. 4-5 
we show that this weak line emission from r > r su b can 
have a dramatic effect onS„ at |u| < HWHM, and thus 
a measurement of S v can provide a new constraint on 
torus models. 

7. SUMMARY 

In this paper, we have shown that the extrapolation 
of the 7 ”blr — L relation found by RM studies of local 
AGN to the most luminous quasars on the sky imply 
^(j'blr) ~ 100 fj, as, with a weak dependence on 2 . Com¬ 
parable angular sizes have previously been resolved in 
young stellar objects using spectroastrometry (Pontoppi- 
dan et al. 2008, 2011), suggesting that spectroastrometry 
is applicable also to the BLR of quasars. 

Assuming that the BLR emission has significant con¬ 
tributions from ordered motions, we calculated the ex¬ 
pected spectroastrometric signal, using a simple model of 
a rotation-dominated BLR. Our calculation suggests that 
the offset between the photocenter of the red-wing and 
blue-wing photons of the broad Ha in luminous quasars 


is detectable (S/N ss 10) with modest time allocations 
on existing telescopes, and therefore can be used to con¬ 
strain telr hr high-L and high -2 quasars. 

This estimate of telr implies a new method to spa¬ 
tially resolve the kinematics of the BLR and directly 
test whether BLR motions are ordered. It also implies 
a new method to contrain Mbh- Since previous meth¬ 
ods to estimate Mbh are based on either RM or on host 
properties, which work best in low-L low -2 quasars, spec¬ 
troastrometry is novel in its ability to estimate Mbh hr 
luminous quasars during the peak of the quasar epoch. 

With next-generation telescopes, BLR spectroastrom¬ 
etry should be routinely detectable for much fainter 
quasars out to z ~ 6, and can be expanded to include 
multiple broad emission lines, including weak lines, at 
an order of magnitude higher velocity resolution than 
achievable with existing telescopes. This will enable de¬ 
mographic and statistical studies of Mbh and BLR prop¬ 
erties as a function of luminosity and redshift, providing 
a compelling science case for next generation telescopes. 
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Table 1 

Properties of the spectroastrometry candidates shown in Fig. 7. Notes: (1) based on mi from §2 and a Richards et al. (2006) SED; (2) 
EW(Ha) from Stern Sz Laor (2012a), EW(Mg n) from Vanden Berk et al. (2001), EW(Pao:) estimated from fig. 1 in Petrov et al. (2012); 

(3) calculated from the line EW and a Richards et al. (2006) quasar SED normalized by L 1450 ; (4) except in 3C 273, all values are 
estimates; (5) from Petrov et al. (2012); (6) based on FWHM(Mg 11 ) from Shen et al. (2011); (7) based on FWHM(Mg 11 ) from Mathews 
Sz Wampler (1985); (8) based on FWHM(C iv) from Paris et al. (2014); (9) based on FWHM(C iv) from Shen et al. (2011); (10) average 
FWHM of other objects; (11) based on FWHM(Mg 11 ) from De Rosa et al. (2014) 











